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A ballistic pn junction in suspended graphene with split bottom gates 
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We have developed a process to fabricate suspended graphene devices with local bot- 



torn gates, and tested it by realizing electrostatically controlled pn junctions on a 
suspended graphene mono-layer nearly 2 /zm long. Measurements as a function of 
gate voltage, magnetic field, bias, and temperature exhibit characteristic Fabry-Perot 



oscillations in the cavities formed by the pn junction and each of the contacts, with 
transport occurring in the ballistic regime. Our results demonstrate the possibility 



CX2 

to achieve a high degree of control on the local electronic properties of ultra-clean 



suspended graphene layers, a key aspect for the realization of new graphene nanos- 
tructures. 
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Technical developments in device fabrication are essential to perform transport exper- 
iments revealing the intrinsic electronic properties of graphene. Suspended graphene de- 
vices 1 ' 2 and devices with hexagonal boron nitride (hBN) as supporting substrate 3 ' 4 provide 
clear examples. The same is true for double-gated devices, in which graphene is not in 
direct contact with any dielectric material. 5 ~ 7 Indeed, these devices have allowed the obser- 
vation of phenomena such as the fractional quantum Hall effect, 8 ~ n new interaction-induced 
symmetry broken states in bilayers, 5 ' 6 and manifestations of ballistic transport. 12 ~ 14 Even 
more advanced experiments would be possible if double-gating on suspended devices could 
be performed locally. In bilayer graphene, for instance, local double gating would allow 
the study of topological confinement, 15 and the realization of fully electrostatically tunable 
pn junctions, of interest to generate or detect light at continuously tunable frequencies in 
the THz to mid infrared range. As an essential step towards the realization of these new 
structures, here we describe a technique to fabricate high-quality suspended graphene de- 
vices with local bottom gates and apply it to the realization of an electrostatically tunable, 
ballistic pn junction in monolayer graphene. 

The fabrication process is illustrated schematically in Fig. 1. The first step consists 
in preparing the bottom gates in the desired configuration -in the present case, a simple 
single strip- on a doped silicon substrate covered with 300 nm SiC>2, by using conventional 
techniques (electron-beam lithography, Ti/Au evaporation, and lift-off). Next, a 450-nm- 
thick layer of polydimethylglutarimide (PMGI)-based lift-off resist (LOR, MicroChem) is 
spun onto the substrate (Fig. 1(a)). LOR resist is chosen because it is not only compatible 
with all subsequent micro-fabrication processes, but also it can be exposed with an electron 
beam and developed away to suspend graphene at the end of the fabrication process. 16 ~ 18 

As a second step, a graphene flake is transferred onto the LOR layer, and positioned on 
to the bottom gate (Fig. 1(b)). To this end, we adapted a technique developed to fabricate 
graphene/hBN heterostructures. 3 Specifically, graphene is exfoliated using an adhesive tape 
and placed on a different substrate, previously coated with a layer of water soluble polymer 
(a 9 wt.% poly(4-styrenesulfonic acid) solution in water) and a layer of PMMA. 19 After the 
desired flake is identified under an optical microscope, the substrate is immersed in water, 
causing the water-soluble polymer to dissolve and the PMMA to float. The floating PMMA 
is retrieved using a plastic support, which is then mounted onto a micro-manipulator under 
an optical microscope. This enables the graphene flake to be transferred onto the LOR, 
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FIG. 1. (Color online) Schematic summary of the fabrication process (drawings not to scale), (a) 
A target substrate with predefined gate electrodes (10 nm Ti/30 nm Au) is covered with a 450- 
nm-thick LOR layer, (b) A graphene flake on a PMM A support is transferred onto LOR, aligned 
to the bottom electrode, (c-d) Graphene is contacted and the underlying LOR exposed with an 
electron beam, to achieve suspension, (e) Optical microscope image of the device (the dashed lines 
indicate the edges of the graphene flake) ; the local bottom gate is visible under the right electrode 
(the bar is 2 /jm long), (e) Schematic top view denoting the regions 1 and 2 (coupled primarily to 
the two different gate electrodes). 



aligned to the bottom gate with a precision of a few microns (Fig. 1(b)). After securing it 
by heating at 105 °C for 40 minutes, the flake is contacted with Ti/Au electrodes (10/60 nm 
thick) denned by conventional electron-beam lithography, metal evaporation, and lift-off (for 
PMMA on LOR, development and lift-off are done using Xylene, 16 ~ 18 at room temperature 
and ~90 °C, respectively). In the final step, the LOR under the graphene layer is exposed 
with an electron beam (Fig. 1(c)) and developed away to achieve the suspension (Fig. 
1(d)). 16 - 18 

We have applied this technique to suspend a 1.8 /iin long graphene monolayer over a 
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FIG. 2. (Color online) (a) Resistance oscillation measured at T = 0.25 K as a function of V% and 
Vi- Along the dashed white line, m changes while 712 is kept fixed at 2.82 x 10 10 cm -2 . The 
data shown in a panel (b), and in Figs. 3 and 4 are taken along this line. The inset shows the 
plot of R as a function of n\ and n%, for the V\ and Vi region corresponding to the parallelogram 
delimited by the dashed line, (b) R(ni) at fixed, positive ri2 shows oscillations for negative n\, 
i.e., when a pn junction is present, (c) Oscillation period Art2, plotted as a function of \/|ra2|. 
Empty squares show data taken at different values of m; the broken line represents the values of 
Ari2 = 2y/TT\n2\/L2 estimated form a simple particle- in-a-box approximation (L2 ~ 800 nm). 



bottom gate that overlaps with about half of the suspended length (see Fig. 1(e)). By 
screening the potential generated by a voltage applied to the conducting Si substrate, the 
bottom gate defines two regions (1 and 2, see Fig. 1(f)), whose carrier density and type 
can be controlled by applying voltages to the doped silicon substrate (Vi) and to the local 
gate itself (V2). Transport measurements as a function of V\ and V2 were performed in a 
Heliox He 3 system to characterize the device at different magnetic field (B), bias (Vdc), and 
temperature (T). Prior to the measurements, the device was annealed at 4.2 K by passing 



a sufficiently large current through the graphene flake. 

Fig. 2(a) shows the resistance R measured at T = 0.25 K as a function of V\ and V 2 . Four 
quadrants can be identified, roughly corresponding to V\ and V 2 having the same or opposite 
sign. When the sign is the same, no pn junction is present in the device: only either electrons 
or holes are accumulated in regions 1 and 2. A pn junction is present between region 1 and 2 
when V\ and V 2 have opposite sign. The borders of the different quadrants are not parallel to 
the V\ and V 2 axis, because of the cross-talk between two gates: V\ does not only change the 
density in region 1 (jii) but also -to a lesser extent- the density in region 2 (712); similarly, V 2 
also influences the density n\. Although, in general, that the density is not spatially uniform 
in regions 1 and 2 (this is obvious when a pn junction is present, in which case the carrier 
density vanishes at the interface between the two regions), accounting as much as possible 
for the effect of the cross-talk is useful to analyze the data. This can be done by looking at 
the gate and magnetic field dependence of the quantized Hall conductance plateaus in the 
unipolar regime, where n\ ~ n 2 (i.e., when the density non uniformity is less pronounced). 
We find m [10 10 cm" 2 ] = 1.0 x V x [V] + 0.35 x V 2 [V] + 0.5 and n 2 [10 10 cm~ 2 ] = 0.2 x V x 
[V] + 1.4 x V 2 [V] — 0.4 (the constants account for the shift of charge neutrality point from 
^1/2 = V; the proportionality terms between 711/2 and V1/2 are in good agreement with the 
estimated geometrical capacitances). The resistance as a function of n\ and n 2 defined in 
this way is shown in the inset of Fig. 2(a). 

When V\ and V 2 are biased with opposite polarity to create a pn junction, the resistance 
doubles as compared to when no pn junction is present (compare, e.g., the resistance for 
n\ < and n\ > in Fig. 2(b)). This shows that the pn junction gives a large contribution to 
the total device resistance, despite the sizable device length (1.8 /xm in total). In particular, 
the pn junction contribution significantly larger as compared to previously studied pn and 
pnp junctions on SiC>2 substrates. 20-23 

Fig. 2(a) further shows that the resistance also oscillates as a function of V\ and V 2 when 
a pn junction is formed, in a way resembling the behavior of graphene pnp junctions on 
a Si/SiC>2 substrate. 19 ' 24 In that case, the oscillations were shown to originate from Fabry- 
Perot interference of Dirac electrons moving ballistically within the (~ 100 nm long) cavity 
defined by the pnp region. 19 ' 24 ' 25 In our device, Fabry- Perot oscillations occur in cavities 
formed by the pn junction and each of the two interfaces with the metal contacts, where 
carriers are also backscattered. 26 ' 27 Our device therefore consists of two Fabry-Perot cavities 
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FIG. 3. (Color online) Magnetic field dependence of R(n{) measured at T = 0.25 K, at fixed n2, 
along the dashed line shown in Fig. 2(a). (a) Plot of R(ni,B), showing characteristic 7r-shift at 
B = B* ~ 20-30 mT. (b) Plot of dG/drix, as a function of n\ and B, shown for comparison with 
similar data reported in the literature. In (a), the broken lines indicate the values of B (B = 0, 
30, 50 mT ) at which the data in (c) are measured. 



connected in series, and the "checkerboard" pattern visible in Fig. 2(a) is a manifestation of 
interference in both cavities. The clear visibility of the oscillations directly in the resistance 
(Figs. 2(a) and 2(b)), without the need of derivating the data, is indicative of the high 
quality of the suspended pn junction. 19 ' 24 

An estimate of the oscillation period An (i.e., the distance in density between two nearest 
resistance peaks or dips) is obtained by imposing that the dynamical phase acquired by an 
electron wave propagating back and forth in the cavity is equal to 2n, i.e. A(2fcpjLj) = 2tt 
(the subscript i — 1, 2 label the region). As kp,i = y/nni, we obtain An, = 2y/irni/Li (note 
that several previous references 19 ' 24 ' 28 reported an incorrect expression, Anj = Ay/Tcrii/Li, 
differing by a factor of 2 from ours). The dotted line in Fig. 2(c) represents the values of 
Ari2 estimated using this formula for region 2 (L2 = 800 nm), and the open squares are 
the experimental values extracted from the most pronounced oscillations measured upon 
changing n2 (a similar result is obtained for region 1). The order of magnitude and the 
trend in the data are well captured by the simple theoretical expression. The experimental 
values, however, are somewhat larger than expected, because the carrier density in the region 
close to the pn junction is lower than the calculated value n^. The lower density causes a 
smaller value of kp, and therefore a smaller phase shift and an additional increase in carrier 
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FIG. 4. (Color online) Energy dependence of the Fabry-Perot interference, (a) Resistance R 
measured at T = 0.25 K as a function of bias Vac an d n\ (at fixed n 2 along a dashed line in 
Fig. 2(a)); the yellow dashed lines are guides to the eye . (b) Two representative curves of R(n\) 
measured at at Vdc = and 5.2 mV, showing a tt shift in the oscillation phase, (c) T-dependence 
of R(ni) measured at Vdc = mV. All data in this figure were taken at B = 40 mT 



density is needed to compensate for this effect. 

The evolution of the oscillation phase upon increasing magnetic field B (see Fig. 3) 
provides further evidence for the Fabry-Perot nature of the interference. 19 ' 24 ' 25 Fig. 3(a) 
shows the 5-dependence of the oscillations upon changing m at fixed n 2 = 2.82 x 10 10 cm -2 
(i.e., by changing V± and V 2 along the dashed line depicted in Fig. 2(a)), which exhibits a it 
phase shift at B = B* rj 20-30 mT (varying n 2 at fixed n\ gives comparable results). Fig. 
3(b) shows the same effect in the derivative of the conductance (G = 1/R) with respect to 
rii, and panel (c) illustrates the occurrence of the phase shift, with three individual slices of 
the color plot shown in (a), taken at B = 0, 30, and 50 mT. 

As discussed for pnp j unctions, 19,24,25 the phase shift originates from the unique prop- 
erties of Dirac electrons, namely the angular dependence of the reflection probability at a 
pn junction, 29,30 and the accumulation of a it Berry phase along momentum-space trajec- 
tories that enclose the origin. 31 ' 32 For a given position of the Fermi energy, the electrons 



contributing predominantly to the Fabry-Perot resistance oscillations are those incident on 
the pn junction with a certain transverse momentum (k y o; the specific value depends on the 
density profile across the junction). 24 ' 25 Upon increasing the perpendicular magnetic field, 
the electron trajectories in the Fabry-Perot cavity are bent, and -in momentum space- 
they eventually enclose the origin. 24 ' 25 When this happens, an additional Berry phase n is 
acquired, causing the phase shift in the resistance oscillations. For ~ 100 nm long pnp 
junctions on substrate, the shift was found to occur at B* ta 2Hk y o/eL m 250-500 mT. 19,24 
Assuming a comparable value of k y o (within a factor of 2-3), this is consistent with our 
observations: the phase shift occurs at an order of magnitude smaller B* ~ 20-30 mT, 
corresponding to an order of magnitude longer cavity. 

Finally, we discuss the characteristic energy scale of the resistance oscillations. Fig. 4(a) 
shows the differential resistance measured as a function of bias Vpc an d density n\. System- 
atically, the position of the resistance peaks shifts linearly upon increasing Vjjc, as expected 
for Fabry-Perot interference. 26,27 The shift is also illustrated by Fig. 4(b), which compares 
measurements taken at Vdc — and 5.2 mV. From both Figs. 4(a) and 4(b), the bias needed 
to shift a maximum of differential resistance into a minimum is approximately 5 meV. We 
have also looked at the energy dependence of the oscillation by changing temperature, and 
found that the oscillations are washed out at about 40 K (f» 3.5 meV). Since, owing to 
the non-uniform charge density, the level spacing in the cavity is somewhat larger than the 
particle-in-a-box value Iivf/ZL ps 2.5 meV (with L ~ 1/im and vf = 10 6 m/s), the energy 
scale found in the experiments is consistent with the simplest theoretical estimate. 

We conclude that the behavior of our device is consistent with the presence of a pn 
junction, and with transport occurring in the ballistic regime over a length comparable to 
the device size (1.8 /im). The measurements therefore confirm that the fabrication technique 
that enables the realization of suspended graphene devices with local bottom gates preserves 
the high quality of the material. In the future, this technology will be applied to realize 
new graphene devices relying on the local control of the electrostatic potential and electric 
field, such as the nano-structures needed for the study of Veselago lensing, 33 of collimation 
of electrons, 34 and of topological confinement. 15 

We gratefully acknowledge A. Ferreira for technical support and financial support from 
the SNF, NCCR MaNEP, and NCCR QSIT. 
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